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Abstract 

The idea is advanced that strong perturbations of an initially equilibrium Bose-condensed 
gas lead to the sequence of nonequilibrium states whose order is inverse to the sequence 
of states arising in the process of the Bose-gas relaxation from an initial nonequilibrium 
state. An approach is described for constructing statistical models of nonequilibrium 
Bose gases. The method is based on the averaging over heterogenous configurations of 
a nonequilibrium system. A statistical model of grain turbulence is suggested. A simple 
model is analyzed consisting of a mixture of two phases, one gauge symmetric and the 
other with broken gauge symmetry. 


1 


1 Introduction 


The description of nonequilibrium systems is notoriously difficult, since such systems are usually 
strongly nonuniform and their dynamical states may quickly vary. However, when it is possible 
to separate in the system dynamics a regime exhibiting specific properties during sufficiently 
long time, longer that the time of fast local oscillations, then it may be admissible to average 
over the local fluctuations and to reduce the consideration to an effective quasi-stationary 
system, whose treatment is essentially simpler than that of the generic nonequilibrium system. 
As examples, we can mention the statistical models of fully developed vortex turbulence [1- 
3], heterophase models of quasi-equilibrium systems [4,5], and effective averaged models of 
nonequilibrium trapped atoms subject to periodic perturbations [6]. 

In the present paper, we describe a general approach for constructing effective quasi- 
stationary models for nonequilibrium systems. The approach is applicable when a nonequi¬ 
librium system exhibits some specific properties during sufficiently long time that is longer 
than the characteristic time of fast oscillations. 

The transfer of a system from an equilibrium state to nonequilibrium states can be done by 
imposing external perturbations. For example, a system can be subject to a time-dependent 
perturbation potential V(t). Then the energy injected into the system by this potential is 
defined as 



In an averaged statistical picture, the injected energy acts similarly to temperature in equilib¬ 
rium systems. Therefore, for a nonequilibrium system of N atoms, it is possible to introduce 
the effective temperature 

rp _ E'inj 

T ‘ ff = teN ' 

which, for brevity, we may denote just as T. 

The outline of the paper is as follows. In Sec. 2, a typical way of creating strongly nonequi¬ 
librium Bose-condensed gases is described and the characteristic parameters of Bose-condensed 
gases, used for creating such strongly nonequilibrium states, are discussed in order to give the 
feeling of the scales involved in experiments with nonequilibrium trapped atoms. In Sec. 3, 
we advance the idea that the procedure of generating a strongly nonequilibrium Bose system 
from an initially equilibrium Bose-Einstein condensate follows the same sequence of states, al¬ 
though in the reverse order, as the process of relaxation of an initially strongly nonequilibrium 
Bose system equilibrating to its condensed state. In Sec. 4, we describe the general method 
of deriving effective quasi-stationary models for nonequilibrium systems. The method is based 
on averaging over heterogeneous configurations arising in the treated nonequilibrium system. 
The consideration is exemplified in Sec. 5 by a statistical model of grain turbulence. In Sec. 6, 
we study a simple statistical model representing a nonequilibrium mixture of two phases with 
different symmetries, one gauge symmetric and the other with broken gauge symmetry. Section 
7 concludes. 

When it does not lead to confusion, we employ the system of units where the Planck and 
Boltzmann constants are set to unity. 
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2 Excitation of trapped Bose-Einstein condensates 

Trapped atoms provide a convenient tool for studying strongly nonequilibrium states of quan¬ 
tum systems. In the process of excitation, an atomic system passes through a sequence of 
qualitatively different states. The system of trapped condensed atoms can be strongly per¬ 
turbed by an external held, when gradually increasing its strength and time of action. 

Trapped Bose atoms in equilibrium at low temperatures form Bose-Einstein condensate in 
the ground state. The condensate cloud in a trap enjoys approximately Thomas-Fcrmi shape, 
with well known properties described in the books [7-9] and reviews [10-19]. In order to break 
the condensate into pieces, it is necessary to impose external perturbations transferring the 
condensate from its ground state to excited states. There are two main ways of imposing such 
external perturbations. 

One possibility is to add to the static trapping potential U( r) an alternating potential 
V(r,t), so that the total trap potential becomes 

U(r,t) = U(r)+V(r,t). (1) 

Another way is to modulate the scattering length a s (t) by means of Feshbach resonance tech¬ 
niques. Both these ways can be used for strongly disturbing condensate [19]. 

Suppose trapped Bose atoms have been cooled down to very low temperatures, when practi¬ 
cally all of them pile down to a Bose-condensed state. And let us apply an external modulating 
perturbation by one of the methods mentioned above. First, at weak perturbation, there ap¬ 
pear elementary collective excitations, that are small deviations from the ground state. Weak 
perturbations also can generate large deviations from the ground state, provided that the mod¬ 
ulation frequency is in resonance with one of the transition frequencies between topological 
coherent modes [20]. The latter are defined as the eigenfunctions of the stationary nonlinear 
Schrodinger (NLS) equation 

- + U(r) + Ar<f> 0 |^ n (r )| 2 <p n (r) = E n ip n (r) , (2) 

2 m 

where N is the number of condensed atoms, assumed to be close to the total number of atoms, 
and 

= 4vr ^ (3) 

m 

is the atomic interaction strength, in which a s is a scattering length assumed to be positive. 
The modes are termed topological, since different modes have different numbers of zeroes, thus, 
topologically different atomic densities. The modes are coherent, being formed by condensed 
atoms characterized by coherent states. 

The known particular example of the topological modes are quantum vortices. If the external 
perturbation rotates the atomic cloud, acting as a spoon, then vortices appear being aligned 
along the imposed axis of rotation. But when the trap modulation does not prescribe a fixed 
rotation axis, then vortices and antivortices arise in pairs or in larger groups [10,20]. The 
explicit experimental demonstration for the appearance of clusters of vortices and antivortices 
was done in Ref. [21]. 

Increasing the strength of the trap modulation generates a variety of coherent modes, need¬ 
ing no resonance conditions because of the power broadening effect [22], Among these numerous 
coherent modes, the basic vortex, with the winding number one, is the most energetically sta¬ 
ble. For a trap with a transverse, ca_i_, and longitudinal, cu 2 , frequencies, the vortex energy can 
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be written [9] as 


COyor 


0.9o;_l 

(FT* 


111(0.81/9). 


where the notation is used for the trap aspect ratio 


Vz 

CJ± 



and for the effective coupling parameter 


g = 4nN -4 


5 


( 4 ) 

( 5 ) 

( 6 ) 


with l± and l z being the transverse and longitudinal oscillator lengths, respectively. Due to the 
large number of atoms N, the effective coupling parameter is large, g > 1. As is seen, the basic 
vortex energy diminishes with the increase of g. At the same time, the transition frequencies 
of other modes, hence their energies, can be shown [20,22] to increase as 

Umn OC (ug) 2/5 (g > 1) . (7) 


This makes the basic vortex the most energetically stable mode. 

When the trap aspect ratio is not too small, the trap can house many vortices, whose 
number can be estimated as 

N„ ~ ^ , (8) 

UJ-vor 

where E in j is the energy injected into the trap by the external pumping. The vortices are 
created due to dynamical instability arising in the moving fluid [23-30]. 

Increasing the strength of the pumping, without imposing a rotation axis, produces a tan¬ 
gle of vortices, which makes the trapped atomic cloud turbulent [31-34]. Increasing further 
either the amplitude of the pumping held or the pumping time leads to the appearance of the 
condensate granulation [35]. 

The energy per particle, injected into the trap by the external perturbation, as is explained 
in the Introduction, plays the role of an effective temperature 





dV (r, t) 

m 


drdt , 


( 9 ) 


where p( r, t) is atomic density. In the case of the periodic in time alternating held V (r, t) ~ 
Acos(cot), the energy per atom, injected during the time interval [fi,f 2 ], takes the form T e ff « 
oj(ti — tf). This makes it possible to represent the crossover lines between different regimes as 
the relation 


A 


T eff 
- 1 2 ) 


( 10 ) 


between the amplitude A of the pumping held and the pumping time. 

The experimental phase diagram on the amplitude-time A — t plane is described in Refs. 
[33-35], where it is shown that with increasing the injected energy, that is proportional to 
the product At, the system passes through the following states: regular superfluid slightly 
perturbed by a weak external held, vortex superfluid with several vortices, turbulent superfluid 
formed by a tangle of many vortices, and granular state with condensate droplets surrounded 
by uncondensed gas. 
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To give the reader the feeling of the typical parameters in the experiments with trapped 
atomic gases, let us mention the corresponding characteristic quantities and experimental data. 

The size of an atomic cloud can be found by solving the nonlinear Shrodinger (NLS) equation 
and calculating the mean-square lengths in the corresponding directions. For a cylindrical trap, 
taking into account that g 3 > 1 , this gives [ 10 , 20 , 22 ] the mean transverse radius 


(2 vg) 1/5 h 
(2tt) 3 / 10 


0.662 (vg) 1/5 l ± 


and the mean axial radius 




{yg) 1/5 lz 

(4 7r ) 3 / 1 ° v /F 


0.468 



( 11 ) 


( 12 ) 


where the standard notations are used for the effective transverse oscillator length l± = 
\Jh/rnojj_ and longitudinal oscillator length l z = a/ K/muj z , and where v is the trap aspect 
ratio. The mean effective cloud radius is 

r 0 = (r 2 ± z 0 ) 1/3 = 0.59 l 0 , (13) 


where the average oscillator length is 


In = 


h 


mujQ 


UJ 0 = (cu 2 ± u z ) 


1/3 


As we see, the actual cloud sizes are noticeably larger than the oscillator lengths because of 
repulsive atomic interactions. These sizes even can be essentially larger than the oscillator 
lengths, when g 3> 1 . 

Knowing the size of the trapped atomic cloud, it is straightforward to find the effective 
cloud volume 

( 7/aW 5 

V eff = nr]_2z 0 = 2wrl = 1.29 -—(14) 
and to estimate the average atomic density in the trap 


P = 


N 


0.775 



(15) 


This shows that, for strongly repulsive atoms, the atomic density p can be much smaller than 
the density N/l g they would have in the absence of repulsive interactions. 

Pair atomic interactions are conveniently characterized by the gas parameter 


7 = p 1,3 a, = 


(16) 


where a = p 1 ^ 3 is mean interatomic distance. The gas parameter 7 is usually small for trapped 
atoms, though can be varied in a wide range by means of the Feshbach resonance techniques. 
Because of the large number of atoms in a trap, the effective coupling parameter g is usually 
large. 

An important quantity, showing whether atoms are in local equilibrium, is the local equili¬ 
bration time ti oc ~ m/hpa s . A perturbed cloud of trapped atoms can, as a whole, be strongly 
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nonequilibrium, while, at the same time, be locally equilibrium. This happens in the situation, 
when the modulation period t mo d = 2 tt/oj of the alternating modulating field, with frequency 
oj, is much longer than the local equilibration time ti oc . 

In experiments [33-35], strongly nonequilibrium states were generated by modulating the 
trapping potential for a trapped cloud of 8 'Rb. The cloud of 87 Rb atoms, of mass m = 1.443 x 
1CT 22 g and scattering length a s = 0.557 x 10 -6 cm, has been cooled down to the temperatures 
much lower than the Bose-Einstein condensation temperature T c = 276 nK, so that the great 
majority of all N = 2x 10 5 atoms have been condensed, the condensate fraction being no = 0.7. 

The trap had cylindrical shape, with the parameters 

u> ± = 1.32 x 10 V 1 , uj z = 1.45 x 10 V 1 , 

l± = 0.74 x 10 _4 cm , l z = 2.25 x 10 _4 cm , 

ujq = 0.63 x 10 3 s _1 , Iq = 1.08 x 10 _4 cm , 

which gives the trap aspect ratio v = 0.11. The effective coupling parameter is g 

The atomic cloud is characterized by the sizes 

r_|_ = 2.27 x 10 _4 cm , z 0 = 1.47 x 10 _3 cm , 

r 0 = 0.42 x 10 _3 cm , V eff = 0.47 x 10“ 9 cm 3 , (18) 

which define the effective atomic density p = 0.43 x 10 15 cm -3 and the mean interatomic 
distance a = 1.32 x 10 ~ 5 cm. The gas parameter is 7 = 0.044. 

The trap potential was modulated by an additional alternating potential V(r,t) [36,37] 
oscillating with frequency u = 1.26 x 10 3 s _1 , which corresponds to the modulation period 
tmod = 0.5 x 10 ~ 2 s. The total modulation time t ext was varied between 0.02 s and 0.1 s. The 
local equilibration time is ti oc = 0.57 x 10 -3 s. Thus, the relations between the characteristic 
times is t ioc <C t mod <C t ext . 

Because of the high atomic density inside the trap, the in situ observation was impossi¬ 
ble. Absorption pictures were taken in the time-of-flight setup, after the times ttof between 
0.015 s and 0.023 s. Restoring the characteristic linear size of granules, corresponding to the 
experimental situation before the free expansion, one gets [35] l g ~ 3 x 10 -5 cm. This is 
in agreement with the theoretical evaluation of the grain size assumed to be of the order of 
the coherence length £ = 1 / \J\.T\pa s . Thus, the relation between the characteristic lengths is 
a s <C a ~ £ ~ l g <C Iq < r 0 . 

The excitation of strongly nonuniform states can also be realized by modulating the scat¬ 
tering length [38,39]. Generally, long modulation times or large exciting amplitudes generate 
the cloud evolution from the appearance of separate vortices to tangled vortex configurations, 
typical of quantum turbulence, and to the granular state. 

To compare the parameters in the experiments with trapped 8, Rb atoms with the typical 
parameters of other experiments with trapped atoms, let us consider the case of 'Li in a 
light trap formed by focused laser beams, as described in Refs. [40,41], The atoms of 'Li 
of mass m = 1.2 x 10 ~ 23 g are prepared, using Feshbach resonance techniques [42], at the 
scattering length a s = 3.2 x 10 ” 8 cm. Lowering down temperature, essentially below the 
critical condensation temperature T c = 200 11 K, almost all N = 3 x 10 5 trapped atoms are 
Bose-condensed, with the condensate fraction no = 0.9. 

The trap characteristics are 

o;_l = 1.48 x 10 3 s “ 1 , l o z = 30.4s _1 , 


(17) 

= 1.96 x 10 4 . 
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l ± = 2.5 x 10" 4 cm , 

l z = 1.7 x 10 3 cm , 


tao = 4.1 x 10 2 s -1 , 

l 0 = 4.7 x 10 _4 cm , 

(19) 

which give the trap aspect ratio v = 0.0207. This means that the trap is quasi-one-dimensional. 
The cloud shape corresponds to the sizes 

r_|_ = 2.62 x 10 4 cm , 

z 0 = 0.88 x 10 2 cm , 


r 0 = 0.84 x 10 _3 cm , 

V eff = 0.38 x 10~ 8 cm 3 , 

(20) 

which define the effective density p = 0.79 x 

10 14 cm -3 and the mean 

interatomic distance 


a = 2.33 x 1CT 5 cm. This density is an order lower than that for 8 'Rb. The lower density 
and the larger spatial size allows to observe the behavior of the atomic cloud in situ. The gas 
parameter is 7 = 1.37 x 10 -3 . The effective coupling parameter is g — 0.48 x 10 3 . 

The effective scattering length 


a s (t) = a B G 


1 


A B 

B(t) - Roo 


can be modulated by varying the magnetic field 


B{t) = B 0 + 5B ■ cos (cut) . 

Here a B Q is a background scattering length far from the resonance field B^, and A B is the 
resonance width. 

In the case, when 5B is much smaller than B 0 , one can write the oscillating scattering length 
in the form 

a s (t ) = a s + 5a s ■ cos (cot) , 

in which 

A B \ _ A BSB 

~ Bo-B^J ’ = aBG (H 0 - B 0 o ) 2 ‘ 

The amplitude of the scattering-length oscillations corresponds to 5a s /a s = 0.2. The fre¬ 
quency u is varied in the range between 157 s ^ 1 and 314 s _1 , which gives the modulation period 
tmod between 2 x 10 ” 2 s and 4 x 10 “ 2 s. The local equilibration time ti oc = 4.5 x 10 ~ 3 s is much 
shorter than the modulation time, ti oc <C t mo d. 

At the beginning, the scattering-length modulation, generates quadrupole mode excitations. 
The amplitude of the oscillations as a function of the applied frequency, allows one to locate 
the resonance curve for the quadrupole mode excitations, as has been shown in Ref. [39]. 

The relation between the characteristic lengths, including the scattering length a s , mean 
interatomic distance a, longitudinal oscillator length l z , and the axial radius of the cloud Zq , is 
such that a s a l z < z 0 . 

The phase diagram on the amplitude-time A — t plane, observed in the experiments with 
a strongly perturbed gas of 8 'Rb, is discussed in Refs. [33-35]. And the detailed results for a 
strongly nonequilibrium gas of 'Li will be published in a separate paper. 


a s = aBG ( 1 
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3 Nonequilibrium condensation versus strong perturba¬ 
tion 

In the process of strong perturbation, trapped Bose gases pass through several nonequilibrium 
states, such as the vortex state with a few vortices, strong vortex turbulence , and granular 
state. Increasing further perturbation should completely destroy the condensate, transferring 
the whole system into a chaotic normal state with no condensate [35]. 

To our understanding, the sequence of dynamic states, appearing in the process of strong 
perturbation of trapped Bose gases, should correspond, although in the reverse order, to the se¬ 
quence of dynamic states, arising in the process of equilibration of an initially strongly nonequi¬ 
librium Bose system to its equilibrium Bose-condensed state. 

The equilibration of weakly interacting Bose systems, from an initial strongly nonequilibrium 
normal state to Bose-condensed state has been studied in a number of publications. Levich 
and Yakhot [43] tried to describe this process by a kinetic equation for the occupation number 
of particles. They found that the time of Bose condensation is infinite, and becomes finite only 
if the presence of germs of the Bose-condensed phase at the beginning of the cooling process is 
assumed. Thus, the single kinetic stage cannot result in Bose-Einstcin condensation. 

Stoof [44,45] employed a functional approach in the frame of the Keldysh formalism for 
deriving a time-dependent Ginzburg-Landau theory and a Fokker-Planck equation for the initial 
stages of nonequilibrium Bose-Einstein condensation. He distinguishes three stages of the 
nonequilibrium condensation. After the Bose gas is quenched into the critical region of the phase 
transition, at the first stage, it can be described by the quantum Boltzmann equation. However, 
such a kinetic equation cannot describe the buildup of coherence in the gas and therefore does 
not lead to a microscopic occupation of the single-particle state. In other words, incoherent 
collisions, governing the Boltzmann kinetic equation, cannot lead to Bose-Einstein condensation 
[43]. To achieve this, a second stage is needed, in which the gas develops the instability toward 
Bose condensation and then coherently populates the ground state by a depletion of the low- 
lying excited states. The assumed instability, developing in the coherent stage, is analogous to 
a dynamic phase transition. After the coherent stage, the gas acquires a highly nonequilibrium 
energy distribution, and equilibrates during the third and final stage. This last stage is again 
of kinetic nature and can be described by the appropriate quantum Boltzman equation for the 
Bogolubov quasiparticles of the Bose-condensed gas. In that way, according to Stoof, there are 
three stages of nonequilibrium condensation: kinetic stage , coherent stage, and relaxation stage 
that also is of kinetic type. 

The kinetic Boltzmann equation was solved numerically by Snoke and Wolfe [46] and by 
Semikoz and Tkachev [47]. In the latter paper, the authors acknowledge the existence of 
three stages, as in the Stoof picture, but treat only two kinetic stages. At the first stage, the 
condensate is absent, but there is a nonzero inflow of particles towards the zero momentum 
state. In the framework of the kinetic equation, there is no condensate at all times if there was 
no condensate initially. Therefore one has to add a seed condensate by hands when switching 
from the first kinetic stage to the final kinetic stage. In this approach, the intermediate coherent 
stage, where the condensate actually emerges by a kind of a phase transition, is omitted. 

The process of Bose condensation has also been considered by numerical solutions of the 
NLS equation assuming the presence of the solution random phase [48] or of external random 
noise [49]. The condensation stages were not clearly distinguished. 

The occurrence of several stages in nonequilibrium Bose condensation of weakly interacting 



Bose gas have been emphasized by Kagan and Svistunov [50,51] and Berloff and Svistunov [52], 
Strongly nonequilibrium Bose gas, after a very short time develops a highly chaotic state, where 
kinetic energy is much larger than the interaction energy. Because of the small nonlinearity, 
the system can be treated as a collection of almost independent modes with random phases. 
The smallness of correlations between the modes implies the absence of any order. This chaotic 
regime of normal (non-superfluid) gas is termed wave turbulence or weak turbulence. Assuming 
that the low-energy modes are macroscopically occupied (so that the occupation numbers are 
much greater than one), it is possible to represent the system by a coherent field with random 
phase. The modes propagate from relatively high to lower energies, and at some time the wave 
turbulence transforms into a regime where short-range coherence starts appearing. This is the 
regime of strong turbulence that cannot be characterized by quasi-independent modes. After 
this intermediate stage, the regime of superfluid turbulence arises, where numerous tangled 
vortices form a random tangle. This could also be called the vortex turbulence. The next stage 
is the process of relaxation of the vortex turbulence to an equilibrium state in a macroscopically 
long time. In this way, one can distinguish four stages: wave turbulence : t strong turbulence , 
vortex turbulence , and relaxation stage. Since from the very beginning, one assumes that the 
low-energy modes are macroscopically occupied, the whole process becomes a crossover, with a 
continuous growth of the low-energy occupation numbers, the lowest of which represents Bose 
condensate. The intermediate regime of strong turbulence is equivalent, in the Stoof picture, 
to his coherent stage. The difference is that Stoof suggests that during this regime a kind of 
phase transition occurs, when the germs of Bose condensate suddenly appear. While in the 
Svistunov et all picture, the low-energy modes are assumed to exist already in the regime of wave 
turbulence, just starting fast growing in the intermediate strong turbulence stage, so that there 
is not a phase transition but a sharp crossover. It has been mentioned [53] that the dynamics 
of the Bose-Einstein condensation is similar to the collapse dynamics of a self-gravitating gas. 

The dynamics of the condensation of a weakly interacting Bose gas in a trap is analogous 
to that of the homogeneous gas [54] and can be characterized by a quantum kinetic equation, 
where the arising condensate comes from the bath of uncondensed atoms [55-57]. 

Zakharov and Nazarenko [58] distinguish four regimes in the dynamics of Bose-Einstein 
condensation. The first is the kinetic stage, when the system is forced into a nonequilibrium 
state corresponding to weak turbulence. The kinetic regime transfers into a strong-turbulence 
state, where the kinetic description breaks down. The second stage is the vortex turbulence, 
where there appear a number of vortices forming a chaotic tangle. In the third stage, the 
system is filled by a well developed condensate with just a few vortices. And the final stage 
corresponds to the relaxation to the equilibrium state. The authors [58] concentrated their 
attention on the strong-condensate regime containing a small number of vortices. They used 
the NLS equation complimented by a term describing forcing and dissipation. 

The condensate dynamics from a nonequilibrium initial state have been studied experimen¬ 
tally [59-62], observing the simultaneous population growth and the development of the phase 
coherence. The regime of vortex turbulence was investigated in experiments [31-35,63] and 
reviewed in Refs. [35,64-67]. 

The regime of vortex turbulence, developing in the process of the nonequilibrium Bose- 
Einstein condensation, is the manifestation of the Kibble-Zurek mechanism [68-70] that has 
been observed [71] in the condensation of trapped 8l Rb atoms. In this picture [68-71], the 
vortex turbulence is preceded by the formation of a nonuniform structure composed of the 
coherent germs, called by Kibble [68] ”cells”, or ”protodomains”, of the condensed phase inside 
the cloud of uncondensed atoms. The order parameters of different cells are random, so that 
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there is no coherence between different protodomains. Such domains move in space and can 
fuse with each other. According to Kibble [ 68 ], it is the cell fusion that creates vortices. 

In order to clearly distinguish these nonequilibrium coherent cells, or protodomains, from 
the static domains occurring in ferromagnets, we shall call these nonequilibrium germs the 
grains. Since their sizes are intermediate between the atomic interaction length and the system 
size, such mesoscopic germs of one phase inside another are analogous to the heterophase fluc¬ 
tuations [4,5]. These coherent grains appear for repulsive atomic interactions, and should not 
be confused with the bound droplet dew forming in a Bose gas with attractive interactions [72], 
Because the shapes and spatial locations of the grains are random, the corresponding strongly 
nonequilibrium phenomenon can be termed grain turbulence. The regime of grain turbulence 
happens between the stages of wave turbulence and vortex turbulence, hence, corresponding 
to the stage of strong turbulence. This is the most difficult regime for theoretical description, 
being strongly nonequilibrium and produced by the chaotic motion of the grains of random 
shapes, which are randomly distributed in space. 

To be more concrete, let us characterize the specific dynamic stages of nonequilibrium 
Bose systems. These stages are connected with the typical length and time scales of statistical 
systems [4,5,12], For dilute gases, the shortest is the interaction length r mt that is much smaller 
than the scattering length a s . Other lengths are: the mean interatomic distance a, correlation 
length £, and the mean free path A, for which we have 


a ~ 


1 



ms ’ 



( 21 ) 


where p is the mean atomic density and 


h 


m 


a/ 47 ipa s 


is sound velocity. For Bose gases, the typical relation between the lengths is 


a s < a < £ < A . 


These lengths are connected with the characteristic velocities: the scattering velocity v s , 
kinetic velocity v a , and the sound velocity s, so that 


h 


ma , 


h 


ma 


h_ 

mt, ’ 


( 22 ) 


with the typical relation 


s < v n < v. 


The kinetic velocity shows the typical velocity of atoms between collisions, while the scattering 
velocity corresponds to the velocity of atoms in the process of their collisions. 

The characteristic time scales, respectively, are as follows. The interaction time 


a s ma s 

Unt ~ ~ —~ 

v q n 


is the time of atomic interactions. The local-equilibration time 


A 

tloc ~ — 


m 


hpa s 


(23) 


(24) 
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is the time during which there develops local equilibrium and there can arise correlated regions 
in space. The heterophase time 


thet ~ tloc ~ ~ —y- (25) 

a s s pa^s 

is the lifetime during which there can exist well correlated regions in space, corresponding to 
protodomaion, or grains. The longer is the relaxation time t re i, after which the system gradually 
relaxes to its equilibrium state during the equilibration time t equ . 

The intervals between the above times define the corresponding dynamic stages through 
which a Bose system passes, being initially prepared in a strongly nonequilibrium state by 
quickly quenching it to the conditions, where the Bose condensate should occur. The first is 
the interaction stage , or dynamic stage, 

0 < t < U n t (interaction stage ) , (26) 

during which atoms interact with each other a few times. At this short initial stage, statistical 
description is not yet applicable. 

The second is the kinetic stage, 

tint < t < ti oc (kinetic stage ) , (27) 

where the system can be described by a kinetic equation. This stage corresponds to the regime of 
wave turbulence, or weak turbulence, when the kinetic energy is much larger than the interaction 
energy, so that the system can be represented as a collection of almost independent modes of 
small amplitudes. The mode independence is clue to their spatial phases being random. Thus, 
the regime of wave turbulence is characterized by three features: large kinetic energy, modes of 
small amplitude, and random spatial phases. There is yet no condensate at this stage. 

The third is the heterophase stage, 


tio C < t < thet (heterophase stage ) , (28) 

when the kinetic description becomes not applicable, since atomic interactions and correlations 
start playing an important role, as a result of which the mode amplitudes grow, and there appear 
mesoscopic well correlated regions of locally condensed atoms. The regions are mesoscopic, 
having the typical sizes of the correlation length £ that is between the scattering length a s 
and the system size. This is the regime of strong turbulence, or grain turbulence. The atoms 
inside each grain are well correlated, forming a condensed droplet, but different grains are not 
necessarily correlated with each other and possess different random phases. The condensed 
grains are surrounded by a gas of normal uncondensed atoms. 

The fourth is the hydrodynamic stage, 


thet. < t < t re i (hydrodynamic stage ) , (29) 

when the mesoscopic condensed grains fuse, forming quantum vortices, according to the Kibble- 
Zurek mechanism [68-71]. A tangle of numerous random vortices arises, creating the regime 
of quantum vortex turbulence. This regime can be described by superfluid kinetic equations 
[73] and superfluid hydrodynamic equations [74,75]. The hydrodynamic stage lasts till the 
relaxation time t re i that depends on the system parameters and geometry [76,77]. 
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After the relaxation time t re i, the vortices decay by mutual recombination and phonon 
emission. For some time, a few vortices survive, but then the system tends to its equilibrium 
Bose-condensed state. This happens at the relaxation stage, 

t re i <t< t equ (■relaxation stage ) , (30) 

where t equ is the equilibration time [76]. 

To illustrate the values of the characteristic parameters discussed above, let us consider the 
experiments [33-35] with trapped 8 'Rb atoms. Then we have the scattering length a s = 0.557 x 
ICR 6 cm, mean interatomic distance a = 1.32 x 1CT 5 cm, correlation length £ 1.823x10 5 cm, 

mean free path A ~ 0.75 x ICR 2 cm, sound velocity s ~ 0.401 cm/s, kinetic velocity v a ~ 0.554 
cm/s, and the scattering velocity v s ~ 13.12 cm/s. This defines the characteristic times: the 
interaction time t int ~ 4.245 x ICR 8 s, local-equilibration time ti oc ~ 0.572 x ICR 3 s, and the 
heterophase time thet ~ 1-87 x ICR 2 s. The interaction time, as is expected, is very short, so 
that atoms quickly pass into the kinetic stage of wave turbulence which lasts around ICR 3 s. 
The duration of the heterophase correlated stage is 1.8 x ICR 2 s, after which the regime of 
vortex turbulence comes into play. 

For other setups, the typical parameters can be rather different. For instance, in the case 
of 'Li, as in experiments [40,41], we have the scattering length a s = 3.2 x ICR 8 cm, mean 
interatomic distance a = 2.33 x ICR 5 cm, correlation length £ ~ 1.775 x ICR 4 cm, mean 
free path A ~ 12.36 cm, sound velocity s ~ 0.495 cm/s, kinetic velocity v a ~ 3.772 cm/s, 
and the scattering velocity v s ~ 2.746 x 10 3 cm/s. From here, we find the interaction time 
ti n t ~ 1.165 x HR 11 s, local-equilibrium time ti oc ~ 4.501 x ICR 3 s, and the heterophase time 
thet ~ 24.97 s. The interaction time is again very short, and atoms quickly pass to the kinetic 
stage of the wave turbulence. But the heterophase correlated stage is quite long, of order of 
ten seconds. Hence, more time is needed for the development of the vortex turbulence, if any. 

The above description corresponds to the relaxation of a Bose gas that initially is prepared 
in a strongly nonequilibrium normal state, after which it tends to its equilibrium condensed 
state. We suggest that when an initially equilibrium condensed Bose gas is subject to strong 
perturbations, it passes through the same dynamic stages, although in the reverse order, with 
the increasing amount of the energy pumped into the system. Thus, an initially condensed 
gas, being perturbed, first goes into a nonequilibrium state, where just a few vortices can arise. 
The next is the hydrodynamic stage, where the vortex turbulence develops. After this, the 
heterophase correlated stage should occur, where the grain turbulence takes place. Finally, 
when all condensate is destroyed, the regime of wave turbulence has to come. The sequence of 
these stages (except the last one) for a strongly perturbed Bose gas of 8 'Rb atoms has been 
observed in experiments [33-35]. The total destruction of the condensate requires very strong 
perturbations that have not yet been reached in experiments. 

The dynamic transitions between the different stages described above are, of course, not 
absolutely sharp, as would be phase transitions in equilibrium systems, but they are rather 
gradual crossovers. 

4 Averaging over heterogeneous configurations 

Although the dynamic transitions between the nonequilibrium stages are crossovers, but inside 
the corresponding temporal intervals the system enjoys qualitatively well defined features. This 


12 



means that it would be admissible to introduce quasi-stationary states as those of a system 
averaged over the appropriate time interval. This can be done as follows. 

Let us assume that a nonequilibrium system can live quite long time, being supported by 
an external perturbation with an alternating modulation potential. If the local equilibration 
time is much shorter than the modulation period, then the system can be treated as quasi¬ 
equilibrium. The difficulty of theoretical description of such a system is in its nonuniformity, 
with randomly arising spatio-temporal fluctuations. When the external pumping potential 
does not impose a spatial symmetry, then the nonuniformities are randomly distributed at each 
snapshot and do not form any ordered structure. Their positions are also random in repeated 
experiments. The nonuniformities are often mesoscopic in space, such that their typical sizes 
are between the interaction radius and the system length. In addition, they are usually of 
multiscale nature, with random sizes and shapes in a dense manifold. Such nonuniformities are 
termed heterogeneous. 

It is possible to show [78] that there exists a mapping between the states of an atomic cloud, 
subject to an alternating modulation, and the states of an atomic system in a random spatial 
external potential. This analogy allows for the estimation of typical parameters characterizing 
the process of nonequilibrium generation. Assume that such a matter has been created. How 
would it be possible to develop a statistical description of this matter? 

Let us consider a snapshot of a heterogeneous matter consisting of the regions of two types, 
whose typical properties can be characterized by some quantities playing the role of local order 
parameters or order indices [79]. For instance, the role of the order parameters can be played 
by densities or some local atomic configurations. Such regions, with different order parameters, 
are analogous to local thermodynamic phases [80,81]. 

The spatial separation of the nonuniformities in a system can be described by employing 
the Gibbs theory of quasi-equilibrium systems [80,81]. The total system volume can be treated 
as a union 

V = V!|JV 2 (31) 

of two parts corresponding to two different spatial regions separated by an equimolecular sur¬ 
face. The convenience of using the equimolecular surface is that it allows for the additive 
representation of observable quantities. For instance, the system volume and the number of 
particles are written as the sums 

V = V 1 + V 2 , N = N, + N 2 . (32) 


Each subvolume is mathematically characterized by the manifold indicator function 




1, r6V„ 

0, r ^ V u 


(33) 


where v — 1,2 enumerates the local phases, that is, the regions with different properties. 

The quasi-equilibrium ensemble, characterizing the system, under a given spatial distri¬ 
bution of different regions, is the pair {"H,/)(£)} of the space of microstates and a statistical 
operator. The space of microstates is the fiber space 




(34) 


of weighted Hilbert spaces corresponding to the states typical of a given region (phase). The 
statistical operator is normalized, so that 


Tr / p{i) = 1 , 


(35) 
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where the trace is taken over the Hamiltonian degrees of freedom and the functional integration 
is over the manifold indicator functions [4], This functional integration characterizes the random 
distribution of the shapes, sizes, and locations of different regions. 

To correctly define the statistical operator, one has to consider a representative statistical 
ensemble taking into account all constraints uniquely describing the system. In addition to the 
normalization condition, the standard constraint is the definition of the internal energy 

E = Tr J p(Z)H{Z)V£ (36) 

as the average of the energy operator H(i). Other constraints can be represented as the 
statistical averages 

C i = Tr[p(()C i (t)Vt (37) 

of the given constraint operators C t (i). The statistical operator is defined as a minimizer of 
the information functional 


mt)} = ^ m^m^i + a 0 


Tr Vi-l 


+ 


+ P 


Tr / p(i)H(i) Vi-E 


Tr f P(()C,({) Vi 


a 


(38) 


in which the first term is the Shannon information and Ao, P, and A* are the Lagrange multipliers 
guaranteeing the validity of the imposed constraints. This principle of minimal information 
yields the statistical operator 

p(0 = exp{— f3H(i)} , (39) 

with the grand Hamiltonian 

ff(£) = ff (£)-£><?, K), (40) 

i 

where /q = —A {F and T — 1/(3 is effective temperature. The inverse normalization factor 


Z — Tr j ex V {-f3H(i)} Vi 


is the partition function. 

The important point is the introduction of the effective Hamiltonian H by the equality 

J exp {-/3H(i)} Vi = exp (-/3H) . (41) 

Then the partition function takes the simple form 

Z = Tre~ pS , (42) 

containing only the trace over the Hamiltonian degrees of freedom. The effective temperature 
T, generally, can include the standard thermal noise and the energy injected into the system, 
as is explained in the Introduction. 
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The geometric weights 


£,(r) Vi 


w v 



(43) 


define the probabilities of the related phases, by construction, satisfying the normalization 
condition 


wi + w 2 = 1 (0 < w u < 1) . 


(44) 


These probabilities are found from the minimization of the grand potential 


!.! = -TlnTre"' 5 " . 


(45) 


The techniques of functional integration over the manifold indicator functions are thoroughly 
explained in Ref. [4]. Below, we shall employ these techniques omitting, for brevity, interme¬ 
diate calculations. 


5 Statistical model of grain turbulence 


As has been mentioned above, the regime of grain turbulence is one of the most difficult for the¬ 
oretical description, if one wishes to consider the details of its nonequilibrium and nonuniform 
spatio-temporal behavior. However, if one is interested in its average statistical properties, it 
is possible to invoke the approach of the previous section. This is admissible, since the local- 
equilibrium time ti oc is much shorter than the lifetime of this regime thet- Here we advance a 
statistical model that can grasp the main average features of the heterophase granular state. 

The granular state can be treated as a heterophase mixture of Bose-condensed grains im¬ 
mersed into the cloud consisting of normal (non-condensed) atoms. Starting with the local- 
equilibrium Gibbs ensemble [80,81], we average over all random heterogeneous configurations, 
as is sketched above, with the related mathematics thoroughly expounded in reviews [4,5,82], 
The resulting effective Hamiltonian of the mixture becomes the sum of two Hamiltonian replicas 

H = H l ^H 2 . (46) 


The Hamiltonian is defined on the fiber space (34). The term H\ acts on the space of microstates 
%i with broken global gauge symmetry U(l), corresponding to the Bose-condensed phase [83], 
while the term H 2 , on the space of microstates 'H 2 with preserved gauge symmetry. 

The Bose-condensed phase is characterized by a representative ensemble [84-87], with the 
grand Hamiltonian 

Hi — Hi — Ho No - Hi Ni - A . (47) 

Here the energy operator is 


Hi = 7i)i / -0 f (r) 


V 2 

2^ + C/ W 


V>(r) dr + 


?// f /S ^ ^ ^ 

+ -H- / ^^(r)^^(r , )<f>(r — r / )'^(r , )r/>(r) drdr' , (48) 

in which U( r) is an external potential, <f>(r) is the interaction potential, and the Bose field 
operators are shifted, according to Bogolubov [88,89], thus, breaking the gauge symmetry, 

^(r) = ?/(r)+^i(r) . (49) 


15 



In this shift, rj( r) is the condensate wave function and t/h( r ) is the held operator of uncondensed 
atoms, these terms being orthogonal to each other, 


J r]*( r)-0!(r) dr = 0 , (50) 

which excludes the double counting. The Lagrange multipliers fi 0 and Hi guarantee the validity 
of the normalization conditions for the number of condensed atoms 


Nq = w x / |r?(r)| 2 dr 


and for that of uncondensed atoms 


Ni — wi / (v4(r)Vh(r))dr 


The last operator A in Eq. (47) is the Lagrange term preserving the condition 

(^i(r)) = 0 

defining the condensate wave function as an order parameter 

V ( r ) = (V>( r )> • 

The grand Hamiltonian of the normal (gauge-symmetric) phase is 


dd 2 — id 2 — ^ 2^2 ■ 


Here the energy operator is 


H -2 = w 2 / (r) — --h U (r) -0 2 (r) dr+ 


+ ^ j , 02( r )t/4( r/ ) < ^( r — r , )t/ ; 2 (r , )t/ ; 2 (r) drdr' . (56) 

The Lagrange multiplier /i 2 guarantees the validity of the normalization for the number of 
atoms in the normal phase, 

N 2 = w 2 J (?/4(r)V’2(r))dr . (57) 

The total number of atoms in the system is 

N = N 0 + N 1 + N 2 . (58) 


The corresponding atomic fractions 


satisfy the normalization 


_ N 0 _ Ni _N 2 

n ° = W' ni = W " 2 = W 


n 0 + ni + n 2 — 1 . 
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The condition of heterophase quasi-equilibrium [4,5,82] leads to the relation 

Vono + Ami = (no + ni)fj, 2 (61) 


between the Lagrange multipliers. 


The geometric weights of the phases are 
This, with the notation 

defined as the 

minimizers of the grand potential. 

Wi = W, 

w 2 = 1 — w 

(62) 

implies the conditions 



on 

d 2 n 

(63) 

dw 

dw 2 


The first of Eqs. (63), with the use of the notations 


Ki = / ($(r) 


K, = 


+ U (r) 


2m ■ - v-,| V’(r)) dr - fio I |r/(r)| 2 dr - fn / (^J(r)^i(r)) dr , 

(V4( r ) _ ^ + ^( r ) ^(r)) dr - fi2 J (</4( r )^2(r)) dr , 

dq = f (ft (r)ip\r')$(r — r , )'?/)(r / )'0(r)) drdr' , 


$2 = / (^2(r)V , 2(r , ) < ^ > ( r _ r')^ 2 (r')^ 2 (r)) drdr' , 


yields the equation 

<l>2 ^ I\ > — I\\ 

W = ----- 

$1 + $2 

for the weight of the Bose-condensed phase. 

The second of Eqs. (63) is the stability requirement imposing the necessary condition 


(64) 


4i + 4 2 > 0 . (65) 

This tells us that the effective atomic interactions must be repulsive in order that the het¬ 
erophase granular mixture would be stable. 

It is worth emphasizing the necessity of nonzero atomic interactions for the Bose-condensed 
system to be stable. The ideal Bose gas is a pathological object that cannot form a stable 
Bose-Einstein condensate because of thermodynamically anomalous particle fluctuations and 
diverging compressibility [12,18,19,77,90,91], although non-condensed ideal Bose gas can be 
stable [92], But in the ideal Bose gas, there can be neither vortices, nor vortex turbulence, nor 
grain turbulence, whose existence is due to the nonlinearity caused by atomic interactions. 

6 Mesoscopic mixture of two phases 

Dealing with the above expressions describing a mesoscopic mixture is a rather involved problem 
requiring the use of nontrivial approximations. Meanwhile, in order to illustrate how the above 
approach works, we consider a simplified model that can exhibit the mesoscopic coexistence of 
two phases, one with broken gauge symmetry and another with the conserved gauge symmetry. 
The choice of this model is done keeping in mind that Bose-Einstein condensation is necessarily 
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accompanied by the global gauge symmetry breaking [83], while the normal, non-condensed 
phase conserves this symmetry. The other advantage of the chosen model is its simplicity 
allowing for a straightforward demonstration of the applicability of the suggested approach. The 
phase with conserved gauge symmetry is called disordered and that with the broken symmetry, 
ordered. The regions of the ordered phase, surrounded by the disordered phase, are mesoscopic 
in the sense that their average sizes are larger than the interaction radius, but much shorter 
than the system linear size. This model can be considered as a cartoon of a system consisting 
of Bose-condensed droplets surrounded by normal uncondensed gas. 

One should not confuse the considered system with mesoscopic thermal disorder with a 
system in a random potential that induces microscopic quenched disorder [93]. The latter 
system can be equilibrium, while the former is only quasi-equilibrium. 

Let us consider an insulating optical lattice, where atoms can occupy several energy levels 
E n . Atoms are assumed to interact through long-range forces, such as spinor or dipole forces 
[94], We start with the standard Hamiltonian in term of the held operators ^(r) that can be 
expanded over atomic orbitals, 

^( T ) = ^2 c nj ( Pnj{r) , (66) 

nj 

where the index n labels energy levels E n and j enumerates lattice cites. Keeping in mind unity 
filling factor, we impose the no-double-occupancy condition 

^ ^ 1 ' C m jC n j 0 . ( 67) 

n 

Taking account of only two lowest energy levels makes it possible to introduce the pseudospin 
operators S° by means of the transformation 

' l/ l.' — 2 + Sj , ~ 2 “ Sj ’ 

Ci,c 2j = s; - iSJ , <4-c„ = S‘ + iS] . (68) 

Assume that the system consists of the random mixture of two different phases. One phase 
has broken gauge symmetry and the other is gauge symmetric. By accomplishing the averaging 
over the manifold indicator functions [4], we come to the effective Hamiltonian 

// = //, ® W 2 , 

H u = w u NE 0 + Y E ^ - w ^ E S i + w l E ( B a s f s i - ^ s i s j) ’ ( 69 ) 

' '•./ 3 ' /./ 

in which A tJ , , I t] are the matrix elements of the atomic interaction potential [95] and 


E 0 = + E2) , Oo = E2 — Ei . 

This Hamiltonian is invariant under the gauge transformation 


S i 




(70) 


(71) 
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in which a is a real number. The pseudospin average, depending on the considered phase, can 
be calculated in two ways, 


(S*) v = Tr n J„S* (is =1,2), 
with different statistical operators 

„ _ exp (~/3H u ) 

Pis — 


Tr«„exp(-/3Jf„ 


(72) 


(73) 


and over different spaces of microscopic states typical of the given phase. The phase with 
broken gauge symmetry corresponds to the nonzero average 


<^>1^0, 

while the phase with the conserved symmetry gives 

< s ;>2 = o . 

The order parameters are defined as 


(74) 


(75) 


N 


S» = 




(76) 


3 =1 


For the phases with the broken symmetry and unbroken symmetry, one has, respectively, 

Si^O, s 2 = 0 . (77) 

To calculate the order parameter si and the phase probabilities w v , we resort to the mean- 
field approximation, with the use of the notations 


A 


B 


fin 


U I + B ’ 

I + B ’ 

U I + B ’ 


1-1 1 ^ 

III 

1 N 

BE-V5 in 

N ^ 3 

*7 £j 

III 

^—i 

(78) 


Then, the minimization of the grand potential ff yields the equations 

2 u + ojiXi — u 2 x 2 


Wi = 


\J (wiSi) 2 + u 2 = w\ tanh 


Au — (1 — b)s\ 

(1 — b)w i 


2 T 


w 2 = 1 - Wi , 

V (whsi ) 2 + uj 2 


(79) 


in which the effective temperature T is measured in units of / + B and 

oji = (1 — b)u , cu 2 — uj — bw 2 x 2 , 


cu 

Xi = — , 

w 1 


x 2 = tanh (j^r) ■ 


(80) 
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We solved these equations numerically for b <C 1, cu <C 1, and for varying u. At T = 0, 
the order parameter si = 1. It monotonically decreases with rising T up to a phase-transition 
point, where Si becomes zero and the system transfers to a disordered phase with unbroken 
gauge symmetry. The phase transition is of second order for u < 0 and u > 3/2, while of first 
order for 0 < u < 3/2. The behavior of w\ and ,Si as functions of T is shown in Figs. 1 and 
2. This model illustrates that, for some system parameters and for sufficiently high effective 
temperature that includes the injected energy, the system, can, first, become a mixture of 
two phases, one with broken gauge symmetry and the other with unbroken symmetry, and at a 
large effective temperature T the whole system transfers to the disordered phase with unbroken 
symmetry, that is, becomes normal gauge-symmetric phase. 

7 Conclusion 

The idea is advanced that strong perturbations of an initially equilibrium Bose-condensed gas 
lead to the sequence of nonequilibrium states whose order is inverse to the sequence of states 
arising in the process of the Bose-gas relaxation from an initial nonequilibrium non-condensed 
state to its equilibrium Bose-condensed state. 

We have described a general approach for constructing statistical models of nonequilib¬ 
rium Bose gases. The method is based on the averaging over heterogenous configurations of a 
nonequilibrium system. 

A statistical model of grain turbulence is suggested, whose general properties are formulated. 
Numerical calculations for a simple model, in the frame of a mean-held approximation, show 
that the grain turbulence can exist in a region of injected energies. Increasing the amount of 
energy, injected into the system above the threshold, leads to the destruction of the regime of 
grain turbulence, so that the system passes to a normal non-condensed state with the preserved 
gauge symmetry. The gauge-symmetric phase represents the wave turbulence of normal, non- 
condensed gas. The transformation of the gauge-broken phase into the gauge-symmetric phase 
plays the role of the transition from the grain turbulence to wave turbulence. 
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Figure Captions 


Fig. 1 . Probability of the ordered phase as a function of the dimensionless effective 
temperature for different parameters u: (1) u = 0; (2) u = 0.1; (3) u = 0.3; (4) u = 0.4; (5) 
u = 0.51; (6) u = 0.75; (7) u — 1; (8) u — 1.5. The dots mark the points of the transition from 
the broken-symmetry phase with a nonzero order parameter to the gauge-symmetric phase with 
zero order parameter. 


Fig. 2. Order parameter as a function of the dimensionless effective temperature for 
different parameters u: (1) u = 0; (2) u — 0.1; (3) u = 0.3; (4) u = 0.4; (5) u = 0.51; (6) 
u = 0.75; (7) u — 1; (8) u — 1.5. The dots mark the points of the transition from the broken- 
symmetry phase with a nonzero order parameter to the gauge-symmetric phase with zero order 
parameter. 
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Figure 1: Probability of the ordered phase as a function of the dimensionless effective temper¬ 
ature for different parameters u: (1) u — 0; (2) u — 0.1; (3) u = 0.3; (4) u = 0.4; (5) u = 0.51; 
(6) u = 0.75; (7) u — 1; (8) u — 1.5. The dots mark the points of the transition from the 
broken-symmetry phase with a nonzero order parameter to the gauge-symmetric phase with 
zero order parameter. 
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Figure 2: Order parameter as a function of the dimensionless effective temperature for different 
parameters u: (1) u = 0; (2) u — 0.1; (3) u = 0.3; (4) u = 0.4; (5) u = 0.51; (6) u = 0.75; (7) 
u — 1; (8) u — 1.5. The dots mark the points of the transition from the broken-symmetry phase 
with a nonzero order parameter to the gauge-symmetric phase with zero order parameter. 
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